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Abstract 

We study a pair of commuting difference operators arising from the elliptic C^-face model. 
The operators, whose coefficients are expressed in terms of the Jacobi's elliptic theta function, 
act on the space of meromorphic functions on the weight space of the Ci type simple Lie 
algebra. We show that the space of functions spanned by the level one characters of the affinc 
Lie algebra sp(4, C) is invariant under the action of the difference operators. 

1 Introduction 

In Ref.[l], one of the authors constructed an L-operator for Belavin's elliptic quantum R-matrix 
[2] acting on the space of meromorphic functions on the weight space of the A n type simple Lie 
algebra. The traces of the L-operator, the transfer matrices, give rise to a family of commuting 
difference operators with elliptic theta function coefficient. In Ref.[3], they are actually equivalent 
to Ruijsenaars' operators [4], which are elliptic extension of Macdonald's g-difference operators [5]. 
The aim of the present paper is to take a step toward a generalization of the above construction to 
the root systems other than the type A. In this paper, we construct a pair of commuting difference 
operators acting on the space of functions on the Ci type weight space. 

In the construction of Refs.[l] and [3], a relation between Belavin's elliptic quantum R-matrix 
and the face-type solution of the Yang-Baxter equation (YBE) [6], especially the intertwining 
vectors [7, 8], played the central role. For the root systems other than type A, it is known no 
vertex-type R-matrices nor the intertwining vectors. Nevertheless, the face-type solutions of the 
YBE are known for all classical Lie algebras and their vector representations [6]. We will utilize 
this type of solution to introduce the difference operators. We take traces (see the section 5) of 
the fused Boltzmann weights to obtain a pair of difference operators (Theorem 1). 

We also show that the space which is spanned by the level one characters of the affine Lie 
algebra Sp(4, C) is invariant under the action of the difference operators (Theorem 2). 

The plan of this paper is as follows. In section 2, we prepare the notation used in the text and 
state the main results. In section 3, we review the C^-face model [6] in the vector representation, 
which was given by a set of functions called Boltzmann weights. In section 4, we introduce the path 
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space, on which the set of Boltzmann weights act naturally as linear maps and thereby explain the 
notion of so-called fusion procedure (see for example Ref. [3] and references there in) . We also give 
a set of formula for fused Boltzmann weights, which leads to the explicit formula of our difference 
operators (Theoreml,(ii)). In section 5, we prove the commutativity of the difference operators. In 
section 6, we prove a property that the difference operators preserve a three dimensional subspace 
spanned by the level one characters of the affine Lie algebra sp(4, C) [9]. In appendix, we give a 
formula of a similarity transformation of the Boltzmann weights. 

Our result can be seen as a type C generalization of Felder and Varchenko's work [10], where 
they showed that the Ruijsenaars system of difference operators can be recovered from the dynam- 
ical R-matrices, which is nothing but the face-type solution of the YBE. 

On the other hand, a BC n generalization of Macdonald polynomial theory is studied by Koorn- 
windcr [11]. In Ref.[12] van Diejen constructed the corresponding family of g-difference operators 
and he studied its elliptic extension in Ref. [13]. He succeeded in constructing two elliptic commut- 
ing operators, one is of the 1st order and the other is of the n-th order, so that they give rise to 
an elliptic extension of difference quantum Calogero-Moser system of type BCi [12] in n = 2. It is 
likely that our operators can be identified with his system with special choice of parameters. We 
hope to report on this issue in the near future. 

Extending this work by van Diejen, Hikami and Komori rescently obtained a general family 
of n-commuting difference operators with elliptic function coefficients [14, 15]. Besides the step 
parameter of difference oprators and the modulus of elliptic functions, the family contains ten 
arbitrary parameters. Their construcion uses Shibukawa-Ueno's elliptic i?-operator [16] togather 
with the elliptic X-operators [17, 18], the elliptic solution to the reflecion equation, and can be 
regarded as an elliptic generalization of Dunkl type operator approach to those systems, which 
have been extensively used by Cherednik [19] (see Ref. [20] for BC n case). It would be interesting 
if one can find an explicit relationship between their approach and ours. 



2 Notation and results 

Let f) be a fixed Cartan subalgebra of the simple Lie algebra g := sp(4,C) and denote by f)* the 
dual space of f). We realize the root system R for (g, f)) as R := {±(si ± £2), ±2ei, ±2e 2 } C f)*. A 
normalized Killing form ( , ) is given by (£?,£&) — \&jk- We will often identify the space f) and its 
dual t)* via the form ( , ). The fundamental weights are given by vj\ — £1, w<i — £\ + £2. Let Vd 
be the set of weights for the fundamental representation L{wd)- We have 

Vl = {±£l,±£ 2 }, V 2 = {±(£l±£ 2 ),0}. (2.1) 

Note that, in these cases, the multiplicity of the weights are all one. 

Fix an elliptic modulus r in the upper half plane Or > and a generic nonzero complex number 
h. Let [u] denote the Jacobi theta function with elliptic nomc p := e 2mT (Sr > 0) defined by 

00 

[u] :=«p 1/8 sin7™ ]J(l-2p m cos27ra + p 2m )(l-p m ). 

m— 1 

This is an odd function and has the following quasi-periodicity 

[u + m] = (-l) m M, [u + rar] = (-i)"^-™™ 2 ^™™^ ( m g z y ( 2 .2) 
,(i) 



Let d, d' be 1 or 2. Then the C 2 type Boltzmann weights of the type (d, d') are given as follows. 

u I is 



For any square \ ^ J (\,ii,v,k e ()*) of weights, the Boltzmann weight Wdd'i ^ ^ 
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given by as a function of the spectral parameter w G C. See the next section for the explicit formula 
for Wu, which are expressed by the Jacobi theta function. 
They satisfy the condition 



w dd , 



A 



K V 



u I = unless [i — A, v — k E 2KPd, k — A, v — /i E 2KP d ' 



and solve the YBE 



P V 

K 



u - v) W dd . 



v - w \W, 



<dd" 



A /i 
P V 

T] V 

a k 



u - wj W d 'd" 
u - w 1 VKrfrf/ 



fl v 

7] K 

A ji 

rj v 



v — w 



u — v 



(2.3) 



The original Boltzmann weights in Ref.[6] are of the type (1,1) in the above terminology. We 
generalized it by the fusion procedure (see the section 4) for the present purpose. 
For A G f)* and pGV d (d=l, 2), we put 

A p := (X,p). 

Theorem 1 Let M d {u) (u G C, d = 1,2) be the following difference operators acting on the space 
of functions on t)* 



(M d (u)f)(\) := J2 W d2 
pev d 



A A + 2hp 
A A + 2hp 



u) Tg/(A), 



where Tp-f(X) := f(\ + 2hp). 

(i) We We M d {u)M dl (v) = M d ,{v)M d {u) (u,v E C,d,d' = 1,2). 

(ii) Lei us define the following difference operators independent of the spectral parameter u 

qjt±p 



M 2 . = J- f'I^±l_^7^7> . t 2/ '] - 2A " " 2/ 'l t 2A " + 2/ ' t A " + " - ^ A "+" + 2/ 'l 



p=±ei 
q=±e2 



[X p+q + h] 



p q 



[6ft] [2A P ] [2X q ] [X p+q + h} [X p+q ] 



Then we have M x {u) = F(u)M 1 ,M 2 (u) = G(u)(M 2 - H(u)), where 

[u] [u + 2h] 2 [u + Ah] 



F(u) :- 



-3ft] 2 [ft] 2 



G(u) :-- 



[u - ft] [u] 2 [u + ft] [u + 2ft] [u + 3ft] 2 [u + 4ft] 
[-3ft] 4 [ft] 4 



(2.4) 



and H{u) :- 



[u + 6ft] [u - 3ft] [2ft] 
[u] [u + 3ft] [6ft] 
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In section 6, we introduce a space of Weyl group invariant theta functions, which are preserved 
by the actions of the difference operators. For (i e ()*, we introduce the following operators S T /s,S/s 
acting on the functions on h*: 

(S Tfj f)(X) := exp [2m ((A, /?) + r(/3, /3)/2)] /(A + r/3), 
(S f)(\) :=/(A + /3). 

They satisfy Heisenberg's relations 

SfiS'y — S'ySfil S T ^S T ~y — S T -yS T f3i S~yS T ft = ^ ^ ' ^ S T S sy (2.5) 

(7,/J,eff)- 

Let Q v ,P V be the coroot and coweight lattice respectively. Let W c Gi(f)*) denote the Weyl 
group for (q, h). Let T/i 14 ' be a space of TL^-invariant theta functions defined by: 



T/i 117 := I f is a holomorphic function over I)* 



S Ta f = S a f = f (VaeQ v )l 

/M) = /(A) (Vw e W) / ■ 

It is well-known that the space is spanned by the level one characters of the affine Lie algebra 
sp(4, C), and the dimension of this space is three. 

Theorem 2 We have 

M d {Th w )dTh w (d=l,2). 
The corresponding facts in the case of A type are proved in Refs.[21] and [3]. 

3 The Ci 1} -face model 

Fix an integer n > 2. We review the definition of the C^'-face model given in Ref.[6]. We realize 
the root system R of the type C n as 

R := {±{e 3 ± e k ), ±2ei 1 1 < j < k < n, 1 < I < n}, 

where {ej}™ =1 is a basis of a complex vector space denoted by f)* with a bilinear form ( , ) defined 

by 

The vector space fj* can be identified with the dual space of a Cartan subalgebra f) of the simple 
Lie algebra sp (2n, C) . The fundamental weights Wj ■ ( 1 < j < n) are given by zuj — E\ + e-2 + ■ • • + £j ■ 
Let V denote the set of weights that belongs to the vector representation L{vjx) of sp(2n,C). We 
have 

V = {±£i, ±£ 2 , ■ • ■ ,±£n}- 

Note that the multiplicity of the weights in V are all one. 
We shall use the following notation frequently, 

V := 2KP, and p := 2hp (p e V). 
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The Boltzmann weights are given by a set of functions of spectral parameter u € C defined 
for any square [ 'j of elements of I)*. Let us denote the functions by W ( ^ ^ u J . They 

V V Hi j \ Hi V I 



satisfy the condition 
W 



A n 

K V 



u I = unless \i — A, v — /x, n — A, v — n € V. 



For p,q,r,s & V such that p + g = r + s, we will write 



P 

s Q7J q =W 

r 



A A + p 

A + s A + p + <f 



They are explicitly given as follows: 
[c — u] [u + ft] 



P 

P [u]p 

P 



P 

p\u\q 

q 
q 

q 



[c] [ft] ' 

[c - u] [A p _ g - u] 
[c] [A P - 9 ] 

[c - u] [u] [Ap_ 9 + ft] 



[c] [ft] [A 



(P ^ ±9), 
(P 7^ ±9), 



..m MjVhgJ^+c-n] [2A p + 2ft]n r#±p [A P+ r + ft] 

/LU ' [c][A P+ , + ft] [2AJ n^ ± jA 9+r ] 



"P 



P 

-P 



c - u] [2X P + h-u] [u] [2A P + ft + c - u] [2A P + 2ftj -pj- 
[c] [2Ap + ft] W ftTI fil [9lJ 11 



[Ap+ g + ft] 
[ 2A p] 9#±p [Ap+g] 



[c] [2Ap + ft] 

The crossing parameter c in the above formulas are fixed to be 

c := -(n + l)ft. 

Proposition 1 The Boltzmann weights (3.1,3.2,3.3,3.4,3.5) enjoy the following properties. 
Initial condition: 



(i 




•) 









(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3.5) 

(3.6) 
(3.7) 



Inversion relation: 



Crossing symmetry: 



(i 


?7 























\ [c + u] [c - u] [ft + u][ft- u] 
" «J = <W ^2 • ( 3 - 8 ) 









(i 


») 




V 





g(A,<c) 

5(m, f) 



k A 



c — u 



(3.9) 
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where we put 

g(X, fj,) := [2fj, p ] Y[ W+q] {^ = X+p, peV). 

Reflection symmetry: 



qjt±p 



W 



X fj, 

K V 



u i = g(\,K)g(K,v) w ( X k 



(3.10) 



Proof . The equation (3.7) is trivial. The two types of symmetries (3. 9), (3. 10) are easily checked 
by the explicit form. In the case of A = v the equation (3.8) is reduced to the following 

E[A P + A r + h + c — u] [X q + X r + h + c + u] 
ieP [x p + x r + h][x q + x r + h] Ar 

{c~u}[c + u}[2X p ][2X q + 2h} ! 

W[2X p + hf x p 

[c + u] [2A p + h + u] [X p + X q + h + c - m] 

[u][2Ap + ft][A p + Aq+fi] 
[c - m] [2A g + ft - u] [A p + A g + + c + tt] 

[w][2A 9 + fi][Ap + A 9 + fi] ' { ' j 

Here we denote by G>p the following function 

0Ap:= _!^ n %i±S (pen (3,2) 

One can find a proof of the equation (3.11) in Rcf.[6](see (3.5) and Lemma 3). The cases /i = 
A + 2p (p G V) are trivial. The remaining cases are easily checked by using the following three term 
identity: 

[u + x]\u — x] [v + y] [v — y] — [u + y] [u — y][v + x] [v - x] 

= [x + y][x - y][u + v][u - v] (3.13) 

(u, v,x,y E C). □ 

We adopted a slightly different formulas (3. 3), (3. 4) from the original ones (see (A.1),(A.2)) in 
Rcf.[6]. In Appendix, we will give a similarity transformation (A. 3), (A. 4) which transforms our 
Boltzmann weights into the original ones. Thus one has a way to prove the YBE for our Boltzmann 
weights, since such a transformation does not destroy the varidity of the YBE. If we follow this 
track, however, we must specify the arguments of the square roots contained in the expressions of 
the original formulas and the transformation. This way of proof may require a rather complicated 
discussion. In this paper, we will give a proof of the YBE for our Boltzmann weights directly 
without using the similarity transformation. In fact, our proof here goes quite parallel to the proof 
given in Ref.[6]. 



Theorem 3 The Boltzmann weights W 
for d = d' = d" = 1. 



A n 

K V 



(3.1,3.2,3.3,3.4,3.5) solve the YBE (2.3) 



G 



Proof. Set 



X(X, p, v, K,a,p\ u, v):=y2w 



v 



P 
a 



n 

K 



w 



X 
P 



P 

n 



w 



p 

V 



Y (A, p, v, k, a, p | u, v) := ^ W ^ * ^ 



?7 v 
a k 



u + v\W 



X p 

1] v 



(3.15) 



and 



Z(X, p, v, n,cr,p \ u, v) := X(X, p, v, K,a,p \ u, v) — Y(X, p, v, K,a,p \ u, v) 



(3.16) 



Regarding Z(X, p, v, k, cr, p \ u, v) as a function of u, we denote it by Z(u). 
The equations (3.7) and (3.8) implies Z(0) = Z(—v) = 0. Since we have 

Z(X, p, v, K,a,p \ u, v) = — ^-y-^-l Zip, A, p,v, K,a\c — u — v, u) 

g(v, k) 



(3.17) 



by (3.9), this shows Z(c— v) = Z(c) = also. Thus we have found the four zeros at u = 0, —v, c, c—v 
of Z(u). By the exactly same argument in Ref.[6] using the quasi-periodicity property of Z(u), 
(3.17) and the following symmetry (this follows from (3.10)) 

~/ x in 9{\p)g{p,o-)g{a,n) 

Z(X, p, is, k, a, p\ u, v) = — — — — Z(X, p, a, K,v,p\ v, u), 

g{X,p)g{p,v)g{v, k) 

we can reduce the proof of the YBE to the following two special cases: 

Z{X,X+p,X+p + q,X+p + q + ?,X + q + ?,X + ? \u,v) = 0, (3.18) 

where r ^ ±p, ±q,p i= ±q and 

Z(X,X+p, X,X+p, X,X+p\u,v) = 0. (3.19) 

In the case of the equation (3.18), each side of the YBE contains only one term, and they are 
manifestly the same. A proof of the last case (3.19) can be found in the original literature [6]. 
However, since the proof is brief and seems to contain some typographical errors, we will describe 
details of it in the following for readers' convenience. 

We will prove Z(X, X+p, X,X+p, A, X+p \u, v) = 0. Regarding Y(X, X+p, X,X+p, A, X+p \u, v) 
as a function of X p we denote it by f(X p ). It reads as 



/(A p ) =Ga p - 



[u] [v] [w] s-^ [Xg + A p + h + U\ [Xg + X p + h + v] [X q + X p + h + w] 



E 



[A g + Ap + h} 3 



+ G 



_! [u\ [v\ [w] [2A P + h-u] [2Ap + h-v] [2Ap + h-w] 
Xp ~\c[ 3 [2Ap + ft]3 

[u] [v\ [w] [2Ap + h + 5] [2Ap + h-v] [2X P + h-w] 



E 



cyclic 



[c] = 



[2Ap + H} 3 



[u] [v] [w] [2A P + H-u] [2Ap + h + v] [2X P + h + w] 

A " 2^ [ c ]3 



cyclic 



[2A P + H} 3 
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where we put w = c — u — v,u = c — u, v = c — v,w = c — w and the summation ^cyclic ^ s over 
the cyclic permutations of the three variables (u, v,w). From the explicit form, one can see that 
X(X, X + p, A, A + p, A, A +p \u, v) = /(— \ p - h). We will prove f(X p ) = /(— X p — h). 
Now consider a function 

. . _ [z + X p + h + u\[z + X p + h + v][z + Ap + h + w] [0]' [2z + 2h] yr [z + X q + H] 
{Z) '- [z + Xp + hf [h] [2z + h] 11 [z + A 9 ] ' 

One sees that $(z) is a doubly periodic function of the periods 1 and r. Its poles are located at 
z = —X p — h, X q (q e V), — | +uj(lu = 0, 5, ^^j 1 )- The pole at z — —X p — h is of the second order, 
and the others are simple. 

Let fi(Xp) (i = 1, 2, 3, 4) denote the i-th term of the above function f(X p ). Since we have 

[x q + x p + h + u\[x q + x p + h + v\[x q + Xp + h + w] 

Ba*{z)dz - [A g + X p + W 3 Xq 

the relation ^Rcs$(z)dz = implies fi(X p ) = a(X p ) + b(X p ) where we set 

a(X p ) := G XP [U] \ V }1 W] V Res $(z)dz, (3.20) 
[cj ^ «=-«+ w 

:= GAp MMJ w ] R CS $( z )dz. (3.21) 

[c| J z=— A p -ft 

Here the summation J2 U 1S over the half periods o> = 0, |, |, i^p. 
From (2.2) and (3.6), we have for w = 0, §, §, 

Res = 1 [A " + I + " + ^ + I + " + ^ + I + " + " ] e^M, (3.22) 

where we put £(0) = £(±) = 0, £(f ) = CC^r 1 ) = c. Combining (3.20), (3.22) and Lemma 3 in 
Rcf.[6], we can verify 

o(Ap) + / 4 (Ap) - f 2 (-X p -h) = -a(-X p -h) + f 2 (Xp) - U(-X p -h)=0. (3.23) 
Set <j>(u) = 4- log[u], then the residue Res <&{z)dz can be expressed as 

z= — A„ — ft 



^cyclic 

+ m + ^ W(-A p + A,) - <f>(-Xp + X q - fi)} J . 



q^±p 



Since <f>(u) is an odd function, we have from (3.21) and (3.24) 

NMM [u\[v\[w\ [2A p + 2ft][2A p ] 



b(X p ) - b(-X p - K) = - 3 



L c] 3 [0]'[h] 2 [2A p + fr] 2 
x {(j)(2X p ) + <f>{2X p + 2K) - 2c/)(2Xp + h)} . (3.25) 
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On the other hand, using the identity (see (3.13)) 

[2A p + h + u\ [2X p + h-v] [2X p + h-w}- [2X p + h-u\ [2X p + h + v] [2X p + h + w] 
[4A P + 2h\ 

and its cyclic permutations of (u,v,w), we have 

h(X P ) M-X p -h)- 3 ^3 [2Xp + h]i - ( 3 - 26 ) 

Now from (3.25) and (3.26), we have 

b(X P ) + h(X p ) = b(-X p -h) + f 3 (-X p - h), (3.27) 
where we used the following identity (Lemma 4 in [6]) 

d>(u + ft) + <j>{u - ft) - 24>(u) = , 

[u\ z [u — h\[u + ftj 

Combining (3.23) and (3.27) we obtained f(X p ) = f(—X p — ft). □ 

4 Path space and fusion procedure 

In the previous section we introduced the Boltzmann weights W(u) of the type (1,1) and proved 
that they satisfy the YBE. In what follows, we treat only the case of n — 2. To construct commuting 
difference operators, we need the general types of the Boltzmann weights Wdd'(u), which we call 
the fused Boltzmann weights. 

First let us introduce the notion of the path space. Let d — 1,2. For any aeC and A,p e I)* 
such that ft-Ae 2hVd, we introduce a formal symbol 

See (2.1) for the notation V\ and Vi. We define the complex vector space 

pi^uy . = f Cg\(u) ■ m - A e 2hv d 

^ d ' x y : otherwise 
for each u G C, and the space of paths from A to fi of the type (c?i, ... ,dk\u\,... ,Uk) 

mi,--- >^fc-ief)* 

The following set 

{sf («i) ® («2) ® ■ ■ ■ ® < fc _ x («*) I Mi - w-i e 2ftP di (1 < i < fc), Mo = A, Mfe = v } 
of paifts forms a basis of the space (4.1). Set also 

V(zu u d l ® • • • ® <*) A := P« ® • • • ® ^ 
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and 



A€f>* 

In the following, we will construct the linear operators 



W dd ,(u - v) : V(m u d w *,) -> ® 

which satisfy the following YBE (d, d',d" = 1, 2) 

(id VK dd <(w - u)) (W dd »(u - to) id) (id W d - d »(t; - w)) 
= {W d , d „(v - w) id) (id W dd »0 - to)) (W dd ,(u - u) id) 

: <„) - otf). (4.2) 

First we define a linear operator W(m\, m\) : V{m\ m\) — > ^(tu^ tu") by 

W(^,^)^W®e»:= ^ w( ^ £ u - w) eft*) e». 

K£t)* V ' 

Put Wu(u - v) := W{mt,ml), then the YBE (4.2) for d = d' = d" = 1 is nothing but (2.3). 
To construct W d d'{u — v) other than W\\{u — v), we will formulate the fusion procedure. Put 



W (m^ m^ • • ■ 



VJ-i 



ml* ®w\)^V(vo\® m^ 



where 



W (to? , n7?)" +1 := id^'- 1 ) W , tu?) id®**" 
: ^(tn^ 1 • • • 

■* ^(tu] 11 • • • nj"' 



J) 



"3 + 1 



We also put 

W (n^ 1 • • • , w? 1 rof • • • tuj' ) 
i<j<l 

:V(m^ • • • ro^ t^ 1 ®---®m v J) ->V{w\^ ■ ■ ■ TO?' m^ •••®^ fc ), 

where 

VK(zu^ •••®G7^,n7^) [j ' fe+jl := id® (j_1) W(n7^ • • • ro^, tup) id® (i " j) 
: ^(n^ 1 • • • ml 3 - 1 w^ 1 • • • w\ k tu"' G7i J+1 • • • w\ l ) 

v ' 

Vival 1 ■■■®w\ 3 - 1 ®wl 3 (3 ml 1 • • • ml k m\ 3+1 • • • m\ l ). 
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We will realize the space Vivo^) as a subspace of V{w\ <g> w\ ). For this purpose, let us 
introduce the fusion projector ir^ by specializing the parameter in W [w\, w\): 

^ : = T^«- ft , ro?) : P(n7 1 ! - R ® tu?) -» P(n7? ® <- s ). (4.3) 

Lemma 1 Tfte space 7r ro ,u (7>(n7"~ R ® ftas a basis {f x+T (u)} re -p 2 given by 

/^+>) :=[V, + %A + >) ® e^f+> - ft) 

+ [A,_ p + h]e X +\u) ® e*+§+> - ft), (4.4) 

wftere p = ±£i, g = ±£2 , and 

#(«) := £ [2A p + 2ft]e*+>) ® e^u - ft). (4.5) 

Proof. For p, q E V\, q ^ ±p, we have 

*«S ~ ft) ® e£f («)) = ^ g p j e* + >) (g) e£f (« - ft) = 0, 



^j(e^(u-ft)®e^| +5 ( U )) 



= I P I ^ U I e A + >) ® e*+§ + > - ft) + \ q g g e* +? (u) ® e*+§+> - ft) 
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-2ft] 



-3ft] [A p _ 



([A p _ 9 + ft]e*+>) ® e^§ + > ft) + [A g _ p + ft]^ + ?( U ) ® e ^+f + > - ft)) 



and 



7^5 (e* + >-ft)®e£ +p -(u)) 

= E ( r H-p ) e " + >) ® e "+^ - n ) 

rePi \ —r J 
_ [-ft] [A p+9 - ft] [A p _ 9 - ft] 



-3ft] [A 



q i^rkr E + 2fi H + >) ® e A +? (- - ft) 



Here we have used the three-term identity (3.13). □ 

Thus we know the subspace Tr^(P(wi~ h i^TUi)\) is naturally isomorphic to the space V(w^)\- 
In the following, we will identify the image Im(7T TO «) C V(tu" ® Tu"~ h ) with the space Vivj^) via 

Proposition 2 Define the operators Wdd'(u — v) by 

W 2 i(u - v) := W{ml ® w^ h , w\), W 12 (u - v) := W(m?, w\ ® w\- h ) (4.6) 
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and 



W 22 {u - v) := W{w u x ® w^~ h , vo\ ® w\~ h ^ 



We have 



W dd ,(u - v){V{wl ® wl,f x ) C V{wl, ® wffl. 



Proof. From the definition of 7r ro - (4.3) and the YBE (2.3), 

W h2 (u - v)W 2 ^{u -v- h)(TT m v ® id) 
= (id ® tt^j)^ 1 ' 2 ^ - « - ft)T4^ 2 ' 3 (u - w). 

Applying this to the definition of W2i(u — v), we get 

W 2 i{u~v){V{w^®wlf x ) CV{w\®w^Y x . 

By a same argument, we have 

W 2 < 3 (u -v + K)W ia {u - «)(id ® 7T ro «) 
= (7r ro? ® id)VF 2 < 3 (u - «)W 1>2 (u - v + h), 

and 



(4.7) 



(4.8) 



Together with the equations (4. 7), (4. 8) and the definition of W 22 (u — v), we obtain 



W 22 (u - v){V{w u 2 ® ts75)5I) C V{w v 2 ® z^. 



□ 



We denote by Wdd'( u — v ) the restricted operators Wdd'( u — v )\ i p(vo ,i ®vo v ^ and introduce their 
matrix coefficients by the following equation 



By the construction, the operators Wdd' {u — v) clearly satisfies the YBE (4.2) in operator form, 

A fi 



and their coefficients Wdd' 



- v satisfies the YBE (2.3). For p,r G V d and s,q G 



TV (d, d' = 1, 2) such that p + g = r + swe write for brevity (as far as confusion does not arise) 



X+p 



s\u\q=W tUf l A + J x+d+n 



(4.9) 



We calculate the coefficients of the operator W 2 \{u) as example. In what follows, we will often 
omit the dependence of g x (u) G V(vo^) on u (the spectral parameter) for brevity. Let p G "Pi. 
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From the definitions of / A and W21 (4.5,4.6) we have 



W 2 i(u) / A ® e A+p - W 21 (u) / A ® e A 



=W 21 {u) I ]T [2Ar + 2ft] e A+? ® e A +? ® e A+g j 



= E 



1 s,tePi 

\s+t=p— g 



where we denote by V q (X; s, t; u) the following function 



[2A r + 2ft] VKn 



A A + f 
A + q A + g + s 



J 



A + f A 

X + q + s X + p 



■ft . 



If q £ Pi such that g 7^ ±p, then the functions V q (X; s,t;u) vanish except for (s,t) 
(p, —q) or (—q,p), and one can easily show that 



Vg{\p, -q;u) = K?(A; -g,p; u) 

[(A + ?) P+9 + ft] [(\ + q)- q - p + ny 

This equation implies that the vector 

V q (X;p, -q; u) e A +p ® e A +§ +? + F 9 (A; -g,p; u) e A +? ® e A+ ? 

is proportional to / A +~ and its coefficient (the both hands sides of (4.10)) is calcurated as 

[u -h][u + ft] [u + 3ft] [2ft] [A 9 _ p -h-u] [2X q + 2ft] 
[-3ft] 2 [ft] 2 [A,_ p - ft] [A 9+p + ft] 

by using the three term identity (3.13). This function is labeled by (see (4.9)) 



q \u\ p (q^ ±p). 
p-q 

Let us consider the term for q = p. For all s E V\ we have from the three term identity 

V P (X; s, —s; u) [u — ft] [u + ft] [u + 3ft] [u + ft] -r-r [A p+r + 2ft] 
[2(A + p) s + 2ft] = [-3ft] 2 [ft] ~ 11 _ 



(4.10) 



m 11 [A p+r + ft] 



(4.11) 



The right hand side of this equation is independent of s e V\ . Thus we see that the vector 

^T/ p (A; S ,- S ;u)e A :| +3 ®e A :| +3 
seVi 

is proportional to f A +| and its coefficient is equal to the right hand side of (4.11), which is labeled 
by 



p\u\p. 
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Here we write all fused Boltzmann weights (the coefficients of the operator W2i(u)). They are 
obtained by the three term identity (3.13). We assume p, q G V\ satisfy p ^ ±q. The common 
factor [u — ft] [u + ft] [u + 3ft] [— 3ft] _2 [ft] _1 is dropped. 



p + q 
q [¥] q 

p + q 



[u + 2ft] 



n P nn n = M [2\ q + 2h][x p - q -n] 
q p y q q [h] [2\ q \ [x p . q + ny 

« 9 ® p = \xm\T ] > (4 - 13) 







[2ft] [A 9 _ p -H-u] [2X q + 2h] 



P-3 



[2ft] [2A 9 - m] [Ap_, - ft] 
[ft] [2A 9 ] [A p+9 + ft] ' 



Next we give the example of W\2- In this case, the common factor [it] [u+2ft] [w+4ft] [—3ft] 2 [ft] 1 
is dropped. To obtain them, we use only the three-term identity (3.13). 



P 

p + q \u\p+q 
P 



[u + 3ft] 



-L. [u + h] [2X P - 2ft] [A,_ p + 2ft] 

q—p [u]q-p — ' 



P 



[ft] [2Ap] [A 



Q-Pl 



p 

\u\q-p - 

q 



[X p -q -2h-u] [X p+q - ft] 



[2A P ] [A 



g-pl 



P -q^]0 [ 2 ^[Ap- 9 -ft-«] [2A g -2ft] 



Q 



[ft] [A g _ p ] [A 



p+q\ 



P [2ft] [2A p -h-u] [X p+q + 2ft] 
p+q\u]q-p= jr, 



-P 



[2X P ] [A 



q-pl 



Finally we give the example of W22 ■ They are equivalent to the Boltzmann weights associated 
to the vector representation of the type B 2 Lie algebra (see Ref.[6]). We write only two cases as 
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example, which is used to define the difference operator M 2 (u). We will drop the common factor 
G(u) (2.4) here. 



(4.16) 



p + q 

4°-H 



Vs = ±£2 



[2A r + 2ft] [2X S + 2ft] [X r+S - 5h] [\ r+s + 2ft] [u + 6ft] [u - 3ft] 



[2A r ][2A s 



[A r+S + ft] 



[u] [u + 3ft] 



(4.17) 



The formulas (4.15), (4.16) and (4.17) together give the explicit form of Md (Theorem 1 (ii)). 



We explain how to calculate the fused Boltzmann weight [IT] . According to the definition 



of the operator W22(u) and the vector (4.5), the coefficient of W 22 (u)f£ ® fx w ^ n res P ec t to 
fx ® /a * s equal to 



1 



i2A^,£, |2V + M1W ' 2 ' 
In this summation, if r is equal to —p, then 

W21 

So that (4.18) can be rewritten as 
W 21 



A A 

X+p X + r 



u)W 21 



X+p X+r 
X X 



u + h 



(4.18) 



A A 
X+p X-p 



u = 0. 



A A 

A+p A+p 

A 



[2A p + 2H] 



u W 21 



X 



X+p X+p 
X X 



X+p X + q 



u)W 21 



X+p X+q 
X X 



u + h 



By means of (4.12), (4.13) and (4.14), this function is equal to 



[u - ft] [u] [u + ft] [u + 2h] [u + 3ft] [u + 4ft] [2h] 



( 



V 



[ft] 



E 



i-m 3 w 

[2X q - 2h] [X p+q + 2h + u] [X p+q — h — u] [A p _q - ft] 



u + ft] [u + 2h] -q [x p+q - ft] [x p+q + 2h] 



[2ft] [-3ft] [x p+q ][x p+q + h] 



\ 



[Ap+g] [A p + 9 - ft] [Ap_ g + ft] 

To obtain the formula (4.17), we use the following lemma. 



q^±p 



) 
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Lemma 2 For any p e Pi, we have 

[U + H] [U + 2H] yr [\ p+q - h] [X p+q + 2h] 

[fi] y- [2A g - 2gj [A p+<? + 2fi + u] [X p+q -h-u] [X p - q - h] 
<?/±p 

_ [it][it + 3fi] [2A r + 2fi] [2A S + 2h] [\ r+s - 5ft] [X r+S + 2fi] [it + 6ft] [it - 3ft] 

~ [6ft] [-3ft] r ^ [2X r ][2X s ] [X r+s }[X r+s + ft] + [6ft] [u + 3ft] ' ( ' ] 

S = ±£2 

Proof . Let f(X p ) be (the left-hand side) — (the right-hand side) of (4.19), regarded as a function 
of X p . It is doubly periodic function of the periods 1, r. Let us show that it is entire. The apparent 
poles of f(X p ) are located at 

X p = X g , X p = X q ±h(p,q e Pi, p + q + 0), A p = 0(p£ Pi). 

Note that the left-hand side of (4.19) is clearly invariant under X q ^> —X q , and the right-hand side 
is W- invariant. In view of the symmetry, it suffices to check the regularity at A p = X q , X p = X q — ft 
and A p = 0. By the three-term identity (3.13), it is easy to see that the residue of f(X p ) at 
A p = X q — ft vanishes. Manifestly, the point A p = X q and A p = is regular. 

Now we have proved that f(X p ) is independent of X p . We will show f(—X q — 2ft) = 0. This can 
be directly checked by using the identity (3.13) twice, and the proof completes. □ 



5 Commutativity of the difference operators 

This section is devoted to the proof of commutativity of the difference operators(Theoreml (i)). 
For t G Pd + Pd' we will introduce the matrices A t (X\u,v), B t (X\v,u) whose index set is I t := 
{(p,q) EPdXPd> \p + q = t} : 



A t (\\u,v)<*$ :=W d 



Bt{X\v,u)%^:=W d , 2 



X X+p 
X X + r 

X X + q 



U W d >2 



v)W d2 



X+p X+t 

X+r X+t 

X+q X+t 

X + s X+t 



\ X X + s 

With these matrices, we can write down both the left and right hand sides as 

M d {u)M d ,(v) = tvA t (X\u,v)T h M d ,{v)M d {u) = ^ tvB t {X\v,u) T ? . 

tev d +v d , tev d +v d , 

Let us also define the matrix Wt(A|u — v) with the same index set: 

The YBE (2.3) implies 

W t {X\u - v)A t {X\u,v) = B t (X\v, u)W t {X\u - v). 

By the inversion relation (3.8), it can be seen that FF t (A|u — v) is invcrtiblc for generic u,v € C. 
It follows that tr A t (X\u, v) = tiB t (X\v,u) for all u,v G C. Hence we have M d (u)M d >(v) = 
M d ,(v)M d (u) for all u,veC. 
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6 Space of Weyl group invariant theta functions 

This section is devoted to the proof of Theorem 2. Let Q v ,P V be the coroot and coweight lattice 
respectively. Under the identification f) = f)* via the from ( , ), these are given by 

Q v = Z2e 1 ®Z2e 2 , P v = Q v + Z(ei + e 2 ). 

Lemma 3 For aZZ /3 e P v and d = 1,2, we have 

[S Tli , M d {u)] = [Sp,M d (u)] = 0. (6.1) 

Proof. Note that if p, q e Pi (g ^ ±p) and /? e P v then (i p+q G Z. By the quasi-periodicity (2.2), 
we have 

[(A + T/3) p+g - ft] _ o nil3 . ah [Xp +q - ft] [(A + f3) p+q - ft] _ [\ p+q - ft] 



e 



[(X + T/3) p+q ] ' [A p+9 ] ' [(A + ^p+J [X p+q ] 

Using these equations, we have for all p£?i 



s rfi n %s^ t ^w = e 2 ^)+^^)/ 2 ) n ^tl^T'i^ +^+a 

^± P [ A p+«] ,J± P P + r/3) P+9 ] 

,2 w i((A,/3)+r(0,/9)/2+2/3 p ft) J| I^±l__^l y ( A + T/? + ^) 

[Ap+gj 



n %s^w/(A), 



and 



Note that 2/3 p ft = (p, (3) etc. Hence we have [S T p, Mi(w)] = [S^, Mi(w)] = 0. In the same way, we 
can see that the principal part of M 2 commutes with S T p and Sp , using the equations 

[(A + r(3) p+q - ft] = ^ t(2f3n+ah) [X p+q - ft] [(A + /3) p+9 - ft] = [A p+ q - ft] 

[(A + Tf3)p +q + ft] ' [Ap +9 +ft]' [(A + ft)p+q + ft] [Ap+g + ft]' 

Using (2.2) it is easy to see that the function 

[2ft] [2A P + 2ft] [2A g + 2ft] [X p+q - 5ft] [X p+q + 2ft] „.„ =6n ^ 
Cm(A) - M^AJ [2A,] [A p+9 + ft] [A p+9 ] G + 9 * 0) 

satisfies C p , g (A + /3) = C M (A + r/3) = C p . q {X) (V/3 £ P v ). This means that S T/3 , SX/3 £ P v ) 
commute with a multiplication by C P , 9 (A). □ 

Lemma 4 For all 7 G P v , we ftawe 

S T1 Th w c Tft 117 , S' 7 Tft w c Th w . (6.2) 
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Proof . Let / S Th w and 7 E P v . Since the bilinear form ( , ) is VK-invariant, we have (S TJ f)(wX) = 
(S' ru ,-i( 7 )/)(A). Using (2.5), we can write this as (SV 7 <SV(u>- 1 ( 7 )- 7 )/)(A), which is equal to SV 7 /(A) 
in view of w' 1 ^) — 7 E Q y ■ In the same way, we can show that {S 1 f)(wX) = (5 7 /)(A). 

Evidently <SV 7 / and S 7 / are holomorphic. For all a E Q v , using (2.5) and (7, a) E Z, it can be 
seen that the operators S ai S Ta commute with S 7 , S r7 . Hence 5 r7 / or 5 7 / are fixed by S Ta and 
S„. □ 

Here we prove Theorem 2. 
Proof of Theorem 2. Let / be any function in Th w . In view of (6.1), we have S a M d f — S TCt M d f = 
M d f for all a E Q v C P v . It is clear from the explicit form of M d that M d f(wX) = M d f(X) for 
all u> E W. _ 

Let us show that the function M d f is holomorphic on fy*. For \i £ h* and z E C, we denote by 
.D* the line in f)* defined by 

!>*:={Aeff |(A,/i) + * = 0}. 

The coefficients of the difference operators M d have their possible simple poles along D + P v + tP v , 
where we put 

D:= |J D°U |J B\ 

P eR + qev 2 -{o} 

and i?+ is a fixed set of positive roots. 

Next we will show that for any function / in Th w , M d f is regular along D. Let us consider 
the meromorphic function g := (jJ peR+ [X p ]j M d f, which is regular along D° := \J peR+ D p . Since 

M d f is VF-invariant, it is clear that g is VF-anti-invariant. This implies that g has zero along D° 
and hence M d f is regular along D°. 

The holomorphy along U 9 ep 2 -{o} 1S somewhat nontrivial. Let p = ±e\, q = ±62- Clearly, 
Mi/ is regular along D p+q . Let us consider the function M 2 f. It suffices to show that the following 
function is regular along D p+q : 

\K+ q ~ h U T fm 1 [gg [ 2A p + 2h \ \ 2 \ + 2fi ] - 5 ^ [ a p+3 + 2ft ] f(X) 
[X p+q + h] v qJ[ ' + [6fi] [2A P ] [2A,] [X p+q + h] [X p+q ] J[ '■ 

We note that, for any TU-invariant function /, we have {T p T q f — /) I^r^ = 0. In view of this, the 

residue of the above function along D p+q is easily seen to vanish. Thus we have proved that for 

any function / in Th w , the functions M d f(d = 1, 2) are regular along D. 
For (5, 7 E P v , we have, by the definitions of S T p, S 1 and (6.1), 

M d f(X + (3t + 7 ) = e- 2 ™« A '« + ^)/ 2 )S r ^ 7 M d /(A) 

= e- 2m « A -' 3 )+ r (' 3 '«/ 2 )M^ T/3> S 7 /(A). (6.3) 

Since S T pS-yf belongs to Th w by (6.2), M d S T f)Sjf is regular along D. Then (6.3) implies that 
M d f is regular along D + (ir + 7. The proof is completed. □ 
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Appendix: Similarity transformation 



Our Boltzmann weights in section 3 and the original form in Ref.[6] are slightly different. The 
original form of type (3.3) and (3.4) are given as follows: 



q 
q 

q 

p [u]-q = 

-p 



[c- 



[c] [h] 



[V* + h ] I A 



p-q 



h] 



1/2 



[A 



p+q 



[c] [X p+q + h] 



[\-q] 2 

— 1 (GxpGxg) 1 ^ 2 



(p?±q), (A.l) 
ip±q). (A.2) 



All the other Boltzmann weights ((3.1), (3. 2), (3. 5)) are the same as the ones we adopted in the 

A /i 
k v 



section 3. We denote these weights by Wjmo 
from those by the following way. We introduce an ordering on the set V as 

E\ < s 2 < ■ ■ ■ < s n < -e„ -<■■■-< -e 2 < 
For A,/i S f)*, such that fi — A = q e 2hV, we define the function s(A, fj,) by 

S (A, M ) :=l[[Xp-q}- 1/2 W-q}- 1/2 . 

pep 

p<q 



u I . Our Boltzmann weights are obtained 



(A.3) 



The relation between the Boltzmann weights W in the section 3 and the ones in Ref.[6] is as follows: 



W 



A fj, 

K V 



s(A, k)s(k, v) 



Wjmo 









(i 


») 




V 





(A.4) 
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